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Abstract 

A convergence criterion of cluster expansion is presented in the case of an abstract poly- 
mer system with general pair interactions (i.e. not necessarily hard core or repulsive). As 
a concrete example, the low temperature disordered phase of the BEG model with infinite 
range interactions, decaying polynomially as l/r'^+-*' with A > 0, is studied. 



1. Introduction 

The abstract polymer gas is an important tool to study the high temperature/low density or low 
temperature phase of many statistical mechanics models. Generally speaking, the abstract polymer 
model consists of a collection of objects (the polymers) which play the role of the particles of the 
gas. These polymers have a given activity and they interact via a hard core pair potential suitably 
defined. Typically, one wants to show that the pressure of this polymer gas can be written in terms 
of an absolutely convergent series if the activities are taken sufficiently small. 

The first example of such a model appeared in [9] where the polymers were finite non overlapping 
subsets of the cubic lattice Tl^ . The authors proved convergence of the pressure via the method 
of Kirkwood-Salsburg equations. Subsequently, the same system studied in [9] was treated in [18] 
and [5] via cluster expansion methods based on tree graph inequalities. 

In [TU] the most general version of this system was given. There, polymers were simply a collection 
of objects with a given activity and interacting through an hard core pair potential introduced via 
a symmetric and reflexive relation in the polymer space. Polymers belonging to this relation were 
called incompatible, and compatible otherwise. The hard core condition was simply to forbid 
configurations of polymers containing pairs of incompatible polymers. Differently from the cases 
considered previously, in which polymers had a cardinality and a size, the Kotecky-Preiss polymers 
were characterized only by the activity. 

In [6] the convergence condition for the Kotecky-Preiss polymer gas was slightly improved and 
the proof was greatly simplified, being reduced to a simple inductive argument, as it was shown 
very clearly in [12] and [19]. In particular, in [19] it has been observed that the Dobrsushin's proof 
works even for more general abstract polymer gases, in which polymers may interact through a 
repulsive soft-core pair interaction. 

Very recently [8] the Kotecky-Preiss and the Dobrushin conditions for convergence of the abstract 
polymer gas with purely hard core interactions were reobtained via the standard cluster expansion 
methods and a new improved condition was given by exploiting an old tree graph identity valid for 
hard core systems due to O. Penrose [14] . 
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In all these works, the basically hard core character of the interaction seemed to be an essential 
ingredient to control the convergence. Exceptions can be found in [7], [TT]. In [7] a contour model 
with interaction (exponentially decaying al large distances) is proposed. However the model is 
rewritten in term of the usual hard core polymer gas where polymers are objects more complicated 
than the original contours. This philosophy has also been pursued in [11] where a one-dimensional 
contour model with long range interaction is rewritten in term of new objects with hard core pair 
interactions. 

It would be of interest to treat also cases in which polymers interact via more general pair 
interactions, e.g., not necessarily repulsive, not necessarily hard core, not necessarily finite range. 
Such abstract polymer model could be a useful tool in the study of spin systems at low temperature 
interacting via infinite range polynomially decaying potential, see e.g. ^13j . 

In this paper we develop a model of abstract polymers (of the type of [TU]) with interactions 
more general than the hard-core. Our polymers interact through a "short distance" repulsive (not 
necessarily hard core) pair potential which is non zero only on pair of incompatible polymers, 
plus an a pair potential with no definite sign (hence it can be attractive), acting only on pairs 
of compatible polymers. We give a condition convergence for the pressure of this gas by using a 
cluster expansion method similar to the one developed in [8]. However, differently from [8j, we 
could not use here the Penrose identity, since our interaction is not purely hard-core. We rather 
used another well known tree graph identity originally proposed in [3] and further developed in [1]. 

The rest of the paper is organized as follows. In section 2 we introduce the model, notations and 
the main result of the paper. In section 3 we give the proof of our result (theorem 1). Namely, in 
subsection 3.1. we present the tree graph identity and show how it can be used to bound the Ursell 
coefficients of the Mayer series of our polymer model. In subsection 3.2 we give the convergence 
argument based on map iterations developed in [8j. In subsection 3.3 we conclude the proof of our 
main theorem. Finally in section 4, as an example, we use theorem 1 to study the low temperature 
disordered phase of the BEG model with infinite range interactions with polynomial decay of the 
type l/r'^+^ with A > 0. 

2. Polymer gas: notations and results 

2.1. The model. 

Let V denotes the set of polymers (i.e. V is the single particle state space). We will assume here 
that P is a countable set. We associate to each polymer 7 E P a complex number (a positive 
number in physical situations) which is interpreted as the activity of the polymer 7. We will denote 
z = {z^l^g-p. 

Polymers interact through a pair potential. Namely, the energy E oi & configuration 71, ... ,7^ 
of n polymers is given by 



where pair potential 1^(7, 7') is a symmetric function in "P x "P taking values in MU {+00}. Observe 
that we don't make any hypothesis on the sign of V^(7i,7j) so this interaction could be for some 
pairs attractive and for other pairs repulsive. 

Fix now a finite set A C V (the "volume" of the gas). Then the probability to see the configu- 





l<i<j<n 



2 



ration (71, . . . , 7„) G A" is given by 

Pro6(7i, . . . ,7„) = J-z^^z^^ • • • z^^e"^i<'<J<"^^'^''^^'^ 



where the normahzation constant 'B\ is the grand-canonical partition function in the vohime A and 
is given by 

Sa(z) = 1 + J;1 %.^72---%„e-^-<.-^(^-^^) (2.2) 

">1 ' (7lv,7n)CA" 

Note that the configurations 71, . . . , 7„ for there exist a pair 7^, 7j such that l^(7i, 7j) = +00 have 
zero probabihty to occur, i.e. are forbidden. So, fohowing the tradition, if a pair (7, j') G V x V is 
such that 1^(7,7') = +00, we wih denote by 7 9^ 7' and say that 7 and 7' are incompatible. 

Since we are admitting non purely repulsive interaction among polymers, we also need to require 
that the potential energy E is stable in the classical sense. This can be achieved by imposing that 
there exists a function -6(7) > such that 

n 

X; ViJ,,J,)>-Y,B{J^) (2.3) 

l<i<j<n 1=1 

for all n G N and all (71, . . . , 7„) E V"'. Note that in the case in which (71, . . . , 7„) contains pairs 
of incompatible polymers the l.h.s. of (2.3) is -|-cxd so this inequelity is trivially satisfied. 
The stabiltity condition immediately implies that Ha is convergent and 

r 1 " 

Ha < 1 + X ^ J2 ^7^^^^^ ^ «^^P { 2^ ^7^^^^^^ } ^ 1^1 ™ exp{z^e^(^)} 



nl 

n>l 



7CA 



7eA 

7eA 



Actually, (2.3) implies that Ha(-z) is analytic in the whole C'^' (|A| is the cardinality of A). 

As we said in the introduction, the usual choice available in the literature is that V{j, 7') takes 
values in the set {0, +00} for all (7, 7') E "P x "P and V{'y, 7) = +00 for all 7 E "P (purely hard 
core pair potential) but we remark that the purely repulsive case (i.e. < l^(7i,7j) < +00 for all 
pairs) has also been considered in [19] and [20] . However, in view of the possible connections with 
the low temperature phase of spin systems with infinite range interactions, we think that the most 
interesting situation treated in the present paper is the case ^(7i,7j) < 0, i.e. when an attractive 
potential, possibly infinite range, is acting among polymers. 

2.2. Results. 

The pressure of this gas, namely log Sa, can be written as a formal series through a Mayer expansion 
on the Gibbs factor exp{— X]i<j<j<„ ^(7*, 7j)}- Namely, a standard calculations (see e.g. [5]) gives 



00 ^ 

^ogEA{z) = Y— (l)'^{-fi,... ,-fn)Z'r^.--Zy„ (2.4) 

with 



n! 

n=l (7i,...,7„)cA" 



(/>'^(7l,...,7n) 



1 if n = 1 

E n (e-^(T-^^) - 1) ifn>2 ^^'^^ 



I, geG„ {i,j}€E. 
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where G„ is the set all connected graphs with vertex set {1,2,... ,n}. We recall that a graph 
g G Gn is a pair g = (Vg, Eg) where = {1, 2, . . . , n} is the set of vertices of g and Eg C C 
{1,2,..., n}} is the set of edges of g. We also recall that g = {Vg, Eg) is connected if for any A, B 
such that AU B = Vg and An B = there exists {i,j} G Eg such That A n {i,j} and 

Bn{i,j}^9. 

The equation (2.4) makes sense only for those z for which the formal series in the r.h.s. of (2.4) 
converges absolutely. To study absolute convergence, we will consider the positive term series 

|logSA|(p) = |'/'^(7i,---,7n)|P7i •••P7n (2-6) 

n=l ■ (7i,...,7„)cA" 

where now G [0, +oo), for all 7 G P and p = {P'jjjeP- Of course | logE\{z)\ < \ logE\\{p) for z 
in the polydisk {\z^\ < p^}^^-p. 

We further define, for each 70 G 7^, a function Hp (p) directly related to (2.6) (a "pinned" sum 
defined in the whole set V) as follows 

np(/^) = E^ E l'/>^(7o,7i,---,7n)K---P7n (2.7) 

n=0 ■ (7i,72,...,7n)eP" 

Clearly, if we are able to show that H^{p) converges, then |logE!A|(/9) and hence |logE!A(2;)| for 
\z-y\ < p'y also converge, since it is easy to check that 

|logSA|(p)<|A|supp^„n;°(p) (2.8) 

70 eA 

To understand the meaning of the series 11^(2;) just observe that its finite volume version Hp (p), 
namely 

00 -. 

^a(p^ = J2^_ E l'^^(7o,7i,---,7n)K---/^7n (2-9) 

n=0 ■ (71,72,. ..,7n)eA" 

is directly related to logSA(pA)- It is immediate to see that 

nr(p) = ^|iogSA|(p) (2.10) 

The main result of the paper is a convergence criterion for the positive series (2.7). Such criterion 
can be considered as a generalization of the Kotecky-Preiss criterion for polymer system interacting 
through a pair potential which is not purely hard core. The criterion can be stated as the following 
theorem. 

Theorem 1 . Let p. : V ^ [0, 00) : 7 1— > /Lt-y be a non negative valued function and let, for each 
7 G 7^, p-y G [0, 00) such that 

P7e-^(^) < /x^e-^7ep-^(T'^)'^^, V7 G (2.11) 
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where -6(7) is the function defined in (2.3) and 



' |e-^(')'»:7,) _ l| = 1 if 00 Jj 

jy(7j,7j)| otherwise 



(2.12) 



Then the series n^o(/o) [defined in (2.7)] converges and satisfies p^^ Yi^^^{p) < fi^^. 

Remark. Observe that in the usual case U hard-core one obtains from theorem [1] the usual 
Kotecky-Preiss condition. We recall however when polymers interact just through a purely repulsive 
potential, one can do better than (2.11). In particular, for the purely hard core case it has been 
shown in [8] that the condition (2.11) can be considerably improved by taking advantage of the 
Penrose tree identity [H], (see also |15j . |19j [8]) valid in the case of purely hard core interactions. 



3. Proof of theorem 1. 



The strategy of the proof is quite similar to the one used in [8J . In particular we use here the very 
same convergence argument for positive series which has been developed in [8]. On the other hand, 
in the present case we cannot use the Penrose identity in order to bound the Ursell coefficients 
|(/)^(7i, . . . , 7„)|, since the pair potential is not purely hard-core (and also not purely repulsive). 
We will rather make use of another well known "tree graph identity" originally proved in [3] (see 

also milKITKle]). 

3.1. Tree graph inequality for 10^(70,71, . . . ,7n)| 

We state the so called tree graph identity [3] , [T] , [1] by using the notations of ^17j and ^16j . We use 
the short notation I„ = {1, 2, . . . , n}. A graph t = {In, Et) G Gn is called a tree if and only if its 
edge set Er has cardinality equal to n — 1. Let us denote by the set of trees with vertex set In- 
In the following whenever [/ is a finite set, \U\ denotes its cardinality. 

Lemma 2 (Tree graph identity). Let Vij, with l<i<j<nbe n{n— l)/2 real numbers, then 
the following identity holds 

E n (^'''^-i) = E n (-^.,)/d^v(t„-i,x„_i)e-^(^"-*"-) (3.1) 

SeGn {i,j}eEg T&T„ {i,j}€Er 

where: 

• t„_i denote a set on n — 1 interpolating parameters tn~i = {ti, ■ ■ ■ ,tn-i) £ [0, 1]"~^; 

• X„_i denote a set of "increasing" sequences of n — 1 subsets, X„_i = Xi, . . . such 
as Vi, Xi C In, we must have Xi C ^j+i, \Xi\ = i and Xi = {1}. 

• i^(Xn_i, t„_i) is a convex decomposition of the potential, explicitly given by 

K(Xn_i,t„_i) = ti{{i,j})...tn-i{{i,jmj (3.2) 

l<i<j<n 
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where 

ti{{i j}) = I ^ ^-^ ^ ^ "■'^'^ ^ ^ ^'^ '"^'^^ 

^ 1 otherwise 

(a pair j} such that i €z Xi and j ^ Xi or vice versa is said to "cross" Xi). 
The measure 

J d/x,(t„_i, X„_i) [■■■]= ■ ■ ■ / ^ ' • ■ ■ ' ^' ' 



(3.3) 



comp. T 



/las totaZ mass egtia/ to one (i.e. it is a probability measure). In (3.3) "X„_i comp. r " 
means that for all i = 1, 2, . . . , n — 1, Xi contains exactly i — 1 edges of r and hi is the 
numebr of edges in r which "cross" Xi 

We want to use (3.1) to bond \4>^{'yi, . . . ,7n)|- This formula is useful when the pair potential 
is not purely repulsive. However, due to the restriction Vij finite (otherwise the r.h.s. of (3.1) is 
not well defined), one in general can apply (3.1) only if the pair potential is finite and absolutely 
integrable, see [1], which is a quite restrictive condition. In particular this rules out a pair potential 
with hard core at short distances which is precisely one of the situations we would like to treat. 

We show here that it is possible to give meaning to r.h.s. of (3.1) even when some among the 
Vij^s take the value oo (the l.h.s. of (3.1) makes sense even in this case). We define a cut-off^ed pair 
potential 

{H ii^ioo 'jj 

(3.4) 
V{'yi , 7j ) otherwise 

Note that, from stability condition (2.3), for any fixed n £ N and (71, . . . ,7„) G V"', there is Ho 
(which depends on n and (71, . . . , 7„) € V"") such that, for all H > Hq and for all X C {1, 2, . . . , n}, 

VH{l^,l,)>-Y.B{^i) (3.5) 

Indeed, if X C {l,2,...,n} is such that {7j}iex does not contain incompatible pairs, then, by 
definition (3.4) and inequality (2.3), it follows 

{i,j}<ZX {i,j}<lX i£X 

If X is such that {'yi}i£x does contain incompatible pairs, then there is at least an edge {k, s} such 
that 1^^(71, 7j) = H so taking 



we have, whenever H > Hq 



{i,j}CX iex 
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So, taking Hq = m.aKxc{i,2,...,n} the inequalities (3.5) are satisfied for all X C {1,2, ... ,n}. 
Now, for any fixed (71, . . . , 7„) G 7^" 

0^(7i,...,7n)= lim Yl n (e-^-(^"^^) - 1) 
5gG„ {i,i}eEg 

We can now use (3.1) for the finite potential Vn and we get 

^ n (e-^^(^-^^) - 1) = 

= hm V n (-F^^(7,,7,)) /d/x,(t„_i,X„_i)e-^«(^"-i'*"-i) 

r€T„ {j,j}eE^ 

where 

-f^if(Xn-i,t„_i) = ^ ti({?,j})...t„_i({z,j})yff(7j,7j) (3.6) 

l<i<j<n 

Now, for fixed r = (!„, £^t-) € ?n and (71, . . . , 7„) G T'" , let us consider the factor 

whin, • • • ,7n) = n \VH,j{li,lj)\ I dMr(t„-l,X„_i)e-^«(^"-i'*"-i) 
{ij}6Er 

The edges C are naturally partitioned into two disjoint sets and E„\E^ where E^ = 

{{i, j} C In '■ 7i 9^ 7j)}- Thus also the edges of the tree r are partitioned into two disjoint sets 
and Er\E^ where E!^ = E^DE^ . So we have 

i/;^(7i,...,7n)= n l^^(7i>7.)l n l^^(7i>7.)l / d//r(tn-i, X„_i)e-^«(^"-i'*"-i) 

(3.7) 

Now, recalling the definition (3.6), we can write 

-f'^i?(X„_i,t„_i) = iC;7^^_^^^(X„_i,t„_i) + irv'^^(X„_i,t„_i) 

where 

-f^f/(l-e)H(Xn-l'*n-l) = X] *l({^i})---^n-l({^'i})^(l-£)H(7i,7i) 

l<i<j<n 

and 

-^Veif(Xn-l5t„_l) = 

. . . tn-li{i,j})VeH{7i,lj) 

l<i<j<n 

where £ > and 

({l-e)H if7j'^7j 

U(l-e)Hi7i,lj) = i 

^ otherwise 
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and 



eH if 7i ^ 

V ("ft , 7j ) otherwise 



The potential ti{{i,j}) . . . tn-i{{i,j})V£Hi7i,7j) satisfies, for H larger that e ^Hq 

{i,j}<ZX iex 
for all X C {1,2, . . . ,n}. This fact implies (see e.g. [1], [Igj, [TT]) that 

n 
i=l 



(3.8) 



The potential i^Tf/^^ ^^jj^ (X„_i, t„_i) is non negative and can be bounded, for r/ > 0, as follows 

^C/(i_,)H(X„-l,t„„i) > ti{{i,j}) ■ ■ ■tn-l{{i,j})U(^i-e)H{li,lj) = 

= Y ti{{i,j}).--tn-i{{i,j}){l-e)H + r] ti{{i,j})...tn-i{{i,j})- 

{id}CEH {i,j}cE^\EH 

-V J2 h{{i,j})...tn-l{{i,j}) > 

{iJjcE^EH 

> E h{{i,j})...tn.i{{i,j}){l-e)H+ Y h{{i,j})...tn-i{{i,j})v-\Er\E^'\v 

{^J}CEH {t,j}cE^\EH 



So we get 



i^c/(i_,)«(X„_i,t„_i) > ti{{hj])---tn-i{{h3])Vi]-\E,\E'^\r^ (3.9) 

l<j<j<n 



where V[j is the positive [H, t], e dependent) pair potential given by 



{l-e)H if{i,i}Gi?f 

if{i,i} Gi?A^^ 
otherwise 



Hence, plugging (3.8) and (3.9) into (3.7) we obtain that w]^(7i, . . . ,7„) can be bounded by 

\Er\E?\ 



n l^(7.,7i)l 

{i,j}(iEr\EH 



l-e 



\E^\ 

n y"dM.(tn-i,x„_i)e-^i<«<.<"*^«*'^'»---*"-«^'^^>)^^^ 



{i,j}<^Er 
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Applying now the tree graph identity (3.1) to the pair potential V^^ one conclude immediately (see 
e.g. ^) that, for all H e [0,+oo) 



^-(l-e)H _ 1 



{i,i}eBr 



(3.10) 



\E?\ 



{i,j}<^EH 



Hence, considering that e^'l^-^^^l |e"'' - l^^^^^l = (e*? - l)I^A£^f I ^e get 



«^^(7i,---,7n)<e+^-i^(^») n 



1 - e 



(e^ - 1) 



^(7., 7,) 



and due to the arbitrarity of i] and e which can be taken as small as we please, and using also that 
^(i-e)H{li-:lj) < ^(7i)7i) foi' ^'^y li 9^ Ij ^-^d any finite H, we obtain, for any H > Hq 

ti^]^(7i,---,7n)<e+S"-i«(-^») n le-^'^^-^^^-ll J] l^(7«>7,)l 

{i,3}^EH {ij}eEr\EH 

which is a bound independent of H. So 



«;-(7i,...,70= lim t«]^(7i, • • • , 7n) < e+^-i J] le-^^^-^^^ - 1 \{ m7i,7i)l 



In conclusion we have that 
where 



(7i,---,7n)|<e+S"-^(^>) ^ n ^(^-^J-) 



^(7i,7i) = < 



f |e-^(T»'Ti) _ i| if7i^7j 
. I^(7j,7j)l otherwise 



and hence also, for n > 1 

-^0 : 



(7o,7i,---,7n)| <e+^"=o^(^») J] J] i^(7i,7i) 



(3.11) 



(3.12) 



where r„ is the set of all trees with vertex set = {0, 1, 2, ... , n}. Inserting (3.12) in (2.7) we get 



"=1 (71,72, •••,7n)e-P" 



reTO {i,j}&E^ 
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< e^?(7o) 



If we pose 



and 



n=l (7i.-,7n)GP" r6T0 {i,j}G£r 



P7 = ^7^ 



5(7) 



(3.13) 



We get 



00 

ni;°(p) = i + E;^E E n i^(7.,7.)P7.---P7. 



(3.14) 



n=l t€TO (7i,...,7„)G-P" {iJjG-Br 



|n|;°(p)<e^(^o)|n|;°(p) 

So the convergence of |n|p (p) implies that of |n|p (/)). 



(3.15) 



3.2. Planar rooted trees and convergence 

We think the trees with vertex set = {0, 1, 2, . . . , n} (i.e the elements of T^) as rooted in 0. We 
define a map m : t m{T) which associate to each labelled tree t £ a unique drawing t = m{T) 
in the plane, called the planar rooted tree associated to r, as follows. 

Given r in T^, place the vertex (the root) at the leftmost position of the drawing. From 
there emerge sq branches ending at the first- generation vertices ii, . . . ,iso- Drawn these vertices 
along a vertical line at the right of the root in such way that the higher has the low label and 
labels increase as we go down along the vertical line (ordering increasing label vertices "from high 
to low"). Then iterate this procedure for the descendants of each first generation vertex (i.e. the 
second generation vertices) ii, . . .is^ and so on... (see figure [T|). 




(a) (b) (c) 

Figure 1: the planar rooted trees associated to the trees (a) with edge set {0, 3}, {1, 3}, {2, 3}, {1, 4}, (b) with 
edge set {0, 2}, {0, 3}, {1, 2}, {2, 4} and (c) with edge set {0, 2}, {0, 4}, {4, 3}, {1, 4}. Observe that (b) and (c) are 
different planar rooted tree 

There is a natural partial order -< among the vertices in a rooted tree. For u,v € t, we say that 
u precedes v and write u < v (or v )^ u) \l the (unique) path from the root to v contains u. If 
{f , u} is an edge of t rooted tree, then either v ~i u 01 u < v. li u < v. u is called the predecessor 
and V is called the descendant. The root has no predecessor and it is the extremum respect to the 
partial order relation ^ in t. For each vertex v of t, we will denote by s„ the branching factor of v 
and we denote by f -"^ . . . , v^" the s^ descendants of f , {v^ being the higher and v^^ being the lower 
in the drawing). 
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Clearly the map r t— >• m(r) = t is many-to-one and the cardinality of the preimage of a planar 
rooted tree t (=number of ways of labelling the n non-root vertices of the tree with n distinct labels 
consistently with the rule "from high to low") is given by 

|{rGrO:m(r)=t}| = ^' , (3.16) 

We denote by = the set of all planar rooted trees with n vertices and by T'^'*'' the set of planar 
rooted trees with maximal generation number k; let also T'^ = Un>oT^ = Ufc>oT''''^ be the set of all 
planar rooted trees. 

Let now fj, : V ^ [0, cxd)^ : 7 i-^ /i^ be a positive valued function defined in V and let, for each 
J £ V, R-y £ [0, 00)^ be defined by the equations 

H-y = R^ip-y{fi), J £V (3-17) 

with 

if^ifl) = 1 + E ^n(7;7l,---,7n)At7i •••/^7„ (3.18) 

for certain functions bn : V^^^ —>■ [0, oo). Denoting Rjipj{fi) = T^o(/") equation (3.17) can be 
visualized in the diagrammatic form 



70 



/^70 ~ ^70 



o 
70 



+ 



70 



-• + 

71 70 




71 



+ 



+ 



72 



•71 
•72 



+ 



•7n 



The sum is over all single-generation rooted trees. In each tree, vertices with open circles with 
subscript 7 represents a factor R^, vertices with bullets with subscript 7 a factor fi-y and vertices 
other than the root must be summed over all possible polymers 7. At each vertex with n descen- 
dants, a "vertex function" bn acts, having as arguments the n + 1-tuple formed by the polymer at 
the vertex and the n polymers associated to the n descendants of that vertex. With this represen- 
tation, the iteration T^(/i) = T(T{fj,)) corresponds to replacing each of the bullets by each one of 
the diagrams of the expansion for T. This leads to planar rooted trees of up to two generations, 
with open circles at first-generation vertices and bullets at second-generation ones. In particular, 
all single-generation trees have only open circles. Notice that the two drawings of Figure [T] appear 
in two different terms of the expansion, and hence should be counted as different diagrams. More 
generally, the A;-th iteration of T involves all possible planar rooted trees up to k generations. In 
each term of the expansion, vertices of generation k are occupied by bullets and all the others by 
open circles. A straightforward inductive argument shows that 

fc-i 



1=0 

where we have denoted R = {i?^}^g-p and 



r4(M) = i27o[j^<)(i?) + $«(i?,/i)] (3.19) 




$W(i?) = V nJ V 6,„(7.;7.i,---,7.-)i?7„i---^7„.J (3-20) 
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Here the product Hui-o ^^^^ vertices of t must be done respecting the partial order of the set 

of vertices in t, i.e. if v >- u the v must be at the right of u in the product. The factor ^^^^^{R, ii) 
has a similar expression but with the activities of the vertex of the k-th generation weighted by /x. 
Here wc agree that 60 (7t)) = 1 and ^0 = 1. We are interested in the ^ ^ 00 limit of (3.20). 

Proposition 2 Let : V ^ [0, co)^ : 7 /i^ he a positive valued function and let, for each 7 G "P, 
-R-y G [0,00)^ be defined by the equations (3.17). Let, V7 eT, Pj & [0, 00) such that pj < Liy. Then 
the series 

*7o(p) := IZn] K{lv;ivi,---,lv^v)p^^, ...p^^,A (3.21) 

is finite for each 70 G P. Furthermore 

P7o^7o(p) < (3-22) 

for each 70 G 7^. 

Proof. By definition (^) = E£o^??(p)- ^y (3.19), the fact that T^^in) = fijo for all A; G N, 
and the assumption p-y < FLy for all 7 G T', we obtain that 

n n 

e=o e=o 

for all n. Thus, since the sequence of partial sums of the series /9^Q$^g(p) is monotonic increasing 
and bounded by p.^^, p^g 

%o{p) converges, and p^o$^o(p) < //^q . □ 
3.3. End of the proof of theorem 1 

We first reorganize the sum over labelled trees appearing in formula (3.14) in terms of the called 

planar rooted trees previously introduced. Namely, recalling that is the set of all planar rooted 
trees with fixed root and n vertices (different from the root), we can rewrite the r.h.s. of (3.14) 
as 

ini;°(p) = i+E;^EE E n ^(^-^.)^7i---P7n (3-23) 

n=l ■ teTO tstO (7i,..,7„)c7'" {i,i}G-B, 

m(r)=t 

Observe that the factor 

E n -^(7i>7j)P7l ■■■P7n 

(71,...,7,0CP" {iJ}eEr 

depends only on the planar rooted tree t = m{T) associated to r (labels of r are dummy indices in 
the sum), i.e. 

E n ^(7i>7i)P7i • • •P7« = n I n E ^(7.,7.op7.. | (3.24) 

(7i,-,7n)CP" {i,j}6£r vhvo ^ i=l%ieV ) 



12 



with the convention that H^i Jl-y ^ep ^ i'lv , Iv^-) P'y^i — 1 when = 0. 

So in conclusion, inserting (3.24) into (3.23) and using also (3.16), we obtain 

nvi~p) = E n Trln E ni.,i.^)~P,^\ (3.25) 

Comparing (3.25) with (3.21) we immediately see that |n|p (p) = (^^^{p) provided 

1 " 

hnir. 71, • • • , In) = -M F{l. 7i) (3.26) 

i=l 

so that 

1 " 

^^{^,) = l + Y^- Yl n ^(^' = e^'^^^ ''^"'"^'^^ (3.27) 

">1 ' (7i,.--,7n)6'P" *=1 

Hence proposition [2] yields the criterion (2.11) for the convergence of the series |n|p (/)) defined 
in (2.7). As a matter of fact, by proposition [21 with the identification (3.26), we have immediately 
that the series |n|!p (/o) defined in (3.14) is finite for each 70 G and p^q\II\^{p) < for each 
7o € V- Now recalling (3.15) and (3.13) we obtain p^(,n^„(p) < n^^. 

4. Example. BEG model with infinite range interactions in the low 
temperature disordered phase 

As an example, we consider the Blume-Emery-Grifhths (BEG) model [2] with infinite range 
interactions in the low temperature disordered phase. The model is defined on the cubic unit 
lattice in d-dimensions Z'^ by supposing that in each vertex x G Z'^ there is a spin variable ax 
taking values in the set {0, —1, +1}. These spins interact via the (formal) Hamiltonian 

H = - E [J-y^^^y + Kccy^rla^y] + D J2 (4.1) 
where Jxy > and K^y G M are summable interactions and we put 

= ^ sup E (J-y + \Kxy\) (4.2) 
^^^^ yez-i 

In the region of parameters 

D> J (4.3) 

the ground state is a = 0. This region is called the disordered phase. If Jxy and Kxy are nearest 
neighbor interactions (or finite range), the low temperature disordered phase can be studied using 
the standard Pirogov-Sinai theory. 

We will make here different assumptions on the interactions Jxy and Kxy Namely, we suppose 
that there exist positive constants c, Ji, A and A' (with < A < A') such that 
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Jxy + \K.y\ < T-^^ V{x, y} G (4.4) 

1-^ y\ 

and 

•^-y ^ 1^ _ l\d+x' \^-y\ ^ \x-l\d+x' (4-5) 

where \x — y\ is the usual nearest neighbor path distance, i.e., \x — y\ is the length of the shortest 
path of nearest neighbors connecting x to y. Due to the assumption (4.5) the low temperature 
phase of the BEG model described by the Hamiltonian (4.1), cannot be studied using the standard 
low temperature Pirogov-Sinai, which explicitly requires finite range interactions. If we further 
assume that the polynomial decay is slow, e.g. by supposing 

A' < 2d + 1 (4.6) 

then this model is not even included in the class of models whose low temperature phase can be 
studied via the extension of the Pirogov-Sinai theory to infinite range interactions given in [l3j . 

We'll show in this section that the partition function of the spin model described by Hamiltonian 
(4.1) can rewritten as the partition function of a polymer system of the type considered in the 
previous sections. Then, using theorem 1, we will prove that, in the disordered phase (4.3) and 
with the assumptions (4.4), (4.5), (4.6), such polymer expansion converges for sufficiently low 
temperatures. 

In order to do that, let us put the system in a finite box A C and let us define, for a fixed 
spin configuration a\ in A, the subset of A given hy P = {x G A : ax ^ 0} . We view this set as 
the union of its connected components, i.e. P = U^^^pi with each set pi C A being connected in 
the sense that for each partition A, B of pi (i.e. Au B = pi and A n i? = 0) there exist x & A 
and y B such that |x — y| = 1. The configuration a\ induces a (non zero) spin configuration 
Sp. on each connected component pi of P which is a function Sp. : pi — > { — 1,+1} : x i— > Sx- 
The pairs Pi = ipi,SpJ are the polymers associated to the configuration a a. By construction the 
correspondence a a {pi, • • • ,Pn} is one to one. The distance between two polymers p = {p, Sp) 
andp= ip,Sp) is the number (i(p,p) = miuxi^p^ y£p\x—y\. Note that if {pi, ... are the polymers 
associated to the configuration a\, then necessarily d{pi,pj) > 2 for all {i,j} C {1, . . . ,n}. 

With these definitions we can rewrite the Hamiltonian of the system in a box AcZ'^ with free 
boundary conditions as (here below P is the inverse temperature) 



n 

pHA{a)= W{p,,pJ) + Y,[pD\p^\-A{pi 



l<j<j<n 1=1 

where 



A{pi) =13 [JxySxSy + Kxy] (4.8) 

{x,y}CPt 

Observe now that to sum over configuration (ta in A is equivalent to sum over polymers configura- 
tions {pi, . . . ,Pn} in A such that n > (n = 0, i.e. no polymers, is the ground state configuration) 
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and d{pi,pj) > 2 for all pairs {i,j} C {1, . . . ,n}. Hence the partition function of the system, at 
inverse temperature (3 and with free boundary conditions, is rewritten as 

Za(/3) =^e-^-^('^^) (4.9) 
= 1 + E;7i E Ppi---/>Pne-^'<^^(^^'^^) (4.10) 



n! 

n>l (pi,---,Pn)ev^ 



where 
and 



Pp = e 



-[l3D\p\-A{p)] ^^j^^^ 



'Pa = {p = {p, Sp) : p C A connected, Sp function from ^?to {— 1,+1}} 
We now extend the definition of W{pi,pj) to all pairs in V as 

' -P E [JxySxSy + K^y] if d{pi,pj) > 2 

Wipi,Pj)=l Hp, (4.12) 

. +00 otherwise 

With these definitions it is immediate to see that r.h.s. of (4.10) can be written as 

Za(/3) = 1 + ^1 Pp,...Pp„e-Si<^<^<"^(P^'P^) (4.13) 

n>l (pi,...,Pn)GPJJ 

which is the partition function of a polymer gas of the type (2.2) in which the polymers are elements 
of the set V defined by 

V = ^p = (jp, Sp) : p (Z IJ^ connected and finite, Sp function from p to {—1, +1}} (4-14) 

with activity given in (4.11) and with incompatibility relation -p p d{p,p) < 2. This pair 
interaction W{pi,pj) is stable in the sense of (2.3). As a matter of fact it is easy to check that, for 
all n e N and all (71, . . . , 7n) G V. 

n 

E w{pi,pj)>-j2B{Pi) 

l<i<j<n 1=1 

with 

B{pi) = pJ\pi\ - A{pi) 

where J is defined in (4.2) and A{pi) is defined in (4.8). Again note that we have to check 
this condition on non intersecting sets of polymers since when (71, . . . ,7^) contains one o more 
incompatible pairs this inequality is trivially satisfied. 

So, by theorem 1, the pressure of this polymer gas (i.e. the free energy of our long range BEG 
model) is absolutely convergent if there exist /Xp such that such that 

g-/3(D-J)|p| <^^g-Epep^'(p,p)Mp^ V7GP (4.15) 
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We choose 



(4.16) 



Hence, inserting (4.16) in (4.15), we obtain that the pressure of such contour gas can be written in 
terms of an absolutely convergent series if, for some a > 



^F{p,p)np < a\p\ 



(4.17) 



By bounding again F(p,p) < 1 whenever p oo p (recall that the short range potential U is in this 
case purely hard core), we get 



pep per 

d(p,p)<l d(p,p)>l 



xe/-" per 



< IpI 2c? sup ^/ip I +m^ ^ \W{p,p)\np 



per 



where d{x,p) = min^gp |a; — y|. Observe now that, by (4.4), |VF(p,p)| < /3Ji|p||p|n ('^'^^^ whenever 
d{p, p) = n. Therefore 

/3e7l 



per 

d(p,x)>l 



n>l per 

d{p,x)=n 



n>l per 



n>l 



pnS„(a:)#0 

where Sn = {y & '^'^ '■ \y\ = n}. An easy calculation show that 



xep 



\Sn\ < 



{2dY 



-n 



So we get 



where 



Hence 



max V \W{p,p)\ Hp < pJ2\p\ sup V \p\ fXp 

per xt^ per 

d{p,x)>l xep 



{2dfJi 



oo ^ 



n=2 



^F(p,p)/Xp< \p\ 

pev 



< \P\ 



hd sup ^ fip] + f3J2 sup ^ \p\ Hp 
2d + /3J2 sup ^ 1^1 Hp < J/3\p\ ^ 1^1 HP 



< 



per 
xep 
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where 

Ji3 = 2d + 13 J2 

Therfore, recalling (4.16), convergence condition (4.17) becomes 

00 

n[e-^(^-^)e"]"2"C„ < a (4.18) 

n=l 

where C„ is the number of connected sets of vertices of with cardinality n containing the origin 
(the factor 2" in l.h.s. of (4.18) counts the number of functions from p to { — 1, +1} when \p\ = n). 
Cn can be easily bounded by C" for some C, e.g. one can take C„ < {Ad)"'. So condition (4.18) 
becomes 

00 

a 

n=l -^/^ 

where 

X = Me-^^^-^^ (4.20) 

Formulas (4.19) and (4.20) imply that 



^n(xe")"<-^ (4.19) 



where 



2u + 1 + V4ii + 1 

For example, taking o: = 1/2 and bounding f{u) < 2u/{2u+l) (we are not looking here for optimal 
estimates), we obtain that convergence occurs if 



8dv^(2d+l + ^J2) 

i.e., for all inverse temperatures /3 > /3o, where Pq is the positive solution of the equation 

/3(D-J) 
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